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Abstract 

We study recently proposed ultraviolet and infrared momentum regulators of the model spaces 
formed by construction of a variational trial wavefunction which uses a complete set of many-body 
basis states based upon three-dimensional harmonic oscillator (HO) functions. These model spaces 
are defined by a truncation of the expansion characterized by a counting number (AT) and by the 
intrinsic scale (hw) of the HO basis; in short by the ordered pair (Af, hw). In this study we choose 
for M the truncation parameter N max related to the maximum number of oscillator quanta, above 
the minimum configuration, kept in the model space. The ultraviolet (uv) momentum cutoff of the 
continuum is readily mapped onto a defined uv cutoff in this finite model space, but there are two 
proposed definitions of the infrared (ir) momentum cutoff inherent in a finite-dimensional HO basis. 
One definition is based upon the lowest momentum difference given by hw itself and the other upon 
the infrared momentum which corresponds to the maximal radial extent used to encompass the 
many-body system in coordinate space. Extending both the uv cutoff to infinity and the ir cutoff 
to zero is prescribed for a converged calculation. We calculate the ground state energy of light 
nuclei with "bare" and "soft" JVJV interactions. By doing so, we investigate the behaviors of the uv 
and ir regulators of model spaces used to describe 2 H, 3 H, 4 He and 6 He with JVJV potentials Idaho 
N 3 LO and JISPI6. We establish practical procedures which utilize these regulators to obtain the 
extrapolated result from sequences of calculations with model spaces characterized by (Af, hw) . 

I. INTRODUCTION 

It has long been suggested that the three-dimensional (3d) harmonic oscillator (HO) provides 
a suitable expansion basis for a straightforward variational calculation of the properties of light 
nuclei. In a traditional variational calculation, a trial wavefunction is selected having a form 
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which aims to exploit all of the important features of the Hamiltonian under investigation, and 
its parameters are adjusted to minimize the energy of the few body system pQ. It is appealing to 
generate a trial wavefunction in a completely systematic manner without regard for the details of 
the Hamiltonian under consideration other than the implementation of exact symmetries. The goal, 
then, is to define a complete set of states for a few-body system and to construct and diagonalize 
the Hamiltonian matrix in a truncated basis of these states. The result of the diagonalization is an 
upper bound to the exact eigenvalue of the complete set. With this method, in contrast to that of 
a pre-chosen trial wavefunction expected to capture the physics, a reliable estimate of the accuracy 
attained can be made with the variational upper bound [2j provided that the trial function is 
constructed using the terms of a systematic expansion set and convergence of the diagonalization 
result (such as a ground-state energy) is observed as the basis is increased. 

The algebra appropriate to generating and using trial wavefunctions, based on 3d HO eigenfunc- 
tions, has been given by Moshinsky [3] and others [4]. The trial functions take the form of a finite 
linear expansion in a set of known functions 

v 

where a{^^ are the parameters to be varied and h v are many-body states based on a summation 
over products of HO functions. The advantage of a HO basis is that it is relatively straightforward 
to construct a complete set of few-body functions of appropriate angular momentum and symmetry; 
examples are given in Refs. [H [5] . The trial function must have a definite symmetry reflecting the 
composition of the bound state: fermions or bosons. This trial function ^>t must be quadratically 
integrable and the expectation value of the Hamiltonian must be finite. The expansion coefficients 
(known as generalized Fourier coefficients in the mathematical literature) depend on the upper 
limit (such as an M defined in terms of total oscillator quanta) and are obtained by minimizing 
the expectation value of the Hamiltonian in this basis. Treating the coefficients a!f^ as variational 
parameters in the Rayleigh quotient [6 ], one performs the variation by diagonalizing the many- 
body Hamiltonian in this basis. This is an eigenvalue problem so the minimum with respect to the 
vector of expansion coefficients always exists and one obtains a bound on the lowest eigenvalue. 
The basis functions can also depend upon a parameter (such as the harmonic oscillator energy hui 
which sets a scale) that then becomes a non-linear variational parameter additional to the linear 
expansion coefficients. Such variational approaches were the standard for calculating properties of 
the trinucleon in the decade following the 1960's [5j [7] and have also been applied to three and 
four body alpha particle models of light nuclei and hypernuclei, see [8HT0]). No-core shell model 
(NCSM) calculations and no-core full configuration (NCFC) calculations with a "bare" potential 
are more recent examples of a variational calculation with a linear trial wave function. Here the 
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basis truncation parameter Af and the HO energy parameter hw are variational parameters |11H13| . 

In such a calculation one would like not only to obtain rapid convergence of the eigenvalue and 
wavefunction but one would like this convergence to be to the exact solution. The functional 
analysis theorems needed for the discussion of the convergence properties of a linear trial function 
are displayed in Appendix A of the article "Variational Techniques in the Nuclear Three-Body 
Problem" [2]. We quote from Section 2.3.5 of the article: "It is shown there [Appendix A] that, 
provided the set of expansion functions is suitably complete [ i.e. complete in the energy norm], 
one will eventually obtain convergence [of the lowest approximate eigenvalue] to the exact value 
[by increasing the basis]. Moreover, if the set is constructed systematically, then in general one 
can expect the convergence to be smooth; indeed, we can often predict the rate at which the 
convergence will occur as we shall discuss... In these circumstances, the numerical convergence of 
the upper bound can provide a useful estimate of the accuracy of the calculation, and one which 
is in practice very much more realistic than that derived from the direct lower bound calculation." 
With a linear trial function, the expectation value W(Af) of any bounded operator W will converge 
provided that the energy converges; and one may estimate the accuracy obtained by watching the 
numerical convergence of W(Af) with increasing Af [H] . Such examinations (albeit for rather small 
basis size compared to those used in this study) are displayed in Ref. |15j . 

The rate of convergence and the number of terms needed for this eventual asymptotic rate to 
"start to behave" is of great practical importance for extrapolation 16] . This question is discussed 
at great length by Delves [2] with general theorems and numerical examples for smooth (e.g. 
attractive Gaussian which is finite everywhere) and non-smooth (e.g. attractive Yukawa which 
has a singularity at the origin) local two-body potentials and a variety of trial functions. As an 
example, Delves derives for the harmonic oscillator basis a convergence rate according to the inverse 
squared power of Af for "nonsmooth" potentials such as a Yukawa; a convergence rate expected to 
be independent of the number of particles. He then demonstrates that the binding energies in the 
truncated expansions of Ref. [5] do follow this power law for both the deuteron of the Reid soft-core 
potential |17j (the archetype of a sum of Yukawas with strong high- momentum components) and 
the deuteron and triton of the (separable) Yamaguchi potential [TH] . This is very slow convergence 
compared to other sets of expansion functions popular in atomic and molecular physics and physical 
chemistry; see Table V of Ref. [2! . Slow convergence hinders progress either because the amount 
of computation needed to reach a desired accuracy is prohibitive, or because too many arithmetic 
operations cause excessive round off error [101 1151 \W\ . Indeed, the slow convergence of systematic 
expansions was likely a contributing factor to the replacement of variational methods by finite 
difference methods (based upon the Faddeev decomposition) in the 1970s to treat the three-nucleon 
bound state problem. 
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In a parallel application of functional analysis to a variational calculation by expansion in a 
basis, specific theorems about the asymptotic rate of convergence for the three-body bound state 
were developed by Schneider for a general basis |20| . The conclusion was that "In any particular 
problem the precise rate will depend on the exact form of the Hamiltonian and the operators 
[which determine the set of basis states] chosen". The practical application of that paper was 
to the hyperspherical harmonics (HH) basis using simple schematic two-body potentials. The 
asymptotic rate of convergence of the three-body binding energy was suggested to converge as the 
inverse fourth power of the maximal grand angular quantum number K for a Yukawa potential 
and exponentially fast in K for a Gaussian potential [20]. These theorems were illustrated by 
explicit HH calculations of the 1970's [3T]. The general expectations of these theorems continue to 
backstop extrapolations in contemporary few-body calculations with modern potentials using this 
HH method [2"2Tf2~T|. As suggested in Ref. [20] , the rate of convergence does not depend on the 
number of particles in the bound state. (This analysis was for 3- and 4-body systems which have 
very high first breakup thresholds-the rate of convergence is, however, expected to depend on the 
first breakup threshold of heavier nuclei). Indeed, contemporary HH analyses of the four-nucleon 
bound state bear out this general expectation, although additional criteria for selecting a reduced 
basis have to be specified, and the authors of [25j demonstrate that the inverse power law in K 
can be higher than four for contemporary "nonsmooth" two- and three-body potentials. 

We are unaware of an application of the theorems proved by Schneider [2D] to the HO basis. 
However, a very up-to-date discussion of the full configuration-interaction (CI) method in a HO 
basis does analyze convergence and gives practical convergence estimates for many-electron sys- 
tems trapped in a harmonic oscillator (a typical model for a quantum dot) [26] . A corresponding 
investigation of light nuclei with another CI method, the NCFC approach [15], provides consistent 
and tested uncertainty estimates for ground state energies. The CI method consists of approxi- 
mating eigenvalues of the many-body Hamiltonian with those obtained by projecting the problem 
onto a finite dimensional subspace of the full Hilbert space and diagonalizing the Hamiltonian in 
this model space [27]. Mathematically, this is analogous to a Ritz-Galcrkin method on the model 
space spanned by the basis functions and the analysis of the energy error is equivalent to analysis 
of the corresponding Raleigh-Ritz calculation sketched earlier. The projection can either take the 
form of an Af defined in terms of total oscillator quanta (called "total-energy-cut space" ) or in the 
single-particle quantum numbers (called "single-particle-cut space" ) . The total-energy-cut space is 
used in this study (see Section II) and the latter single-particle-cut space lends itself more readily 
to approximating a full CI calculation by a coupled cluster approach. The CI approach becomes, 
in principle, exact as Af — > oo with either choice of Af. For this reason the CI approach with HO 
basis functions is sometimes called called "exact diagonalization" . A succinct statement of the 
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equivalence of large-scale diagonalization and the Rayleigh-Ritz variational method can be found 
in the Introduction of Ref. [28] and the full discussion in the monograph Ref . [29] . 




As in the nuclear examples [3J [20] , the asymptotic convergence of the lowest eigenfunction of a 
quantum dot in the "total-energy-cut" model space is slow and it is slow for the analogous reason: 
the singularity of the Coulomb interaction at those points where two or more interparticle distances 
are zero (the Kato cusp condition on the many-body wavefunction [30]). The convergence rate is 
dominated by the singularities in the analytic structure of the solution [31] : it takes many HO 
eigenfunctions to approximate the singularities of the many-body wavefunction due to the two- 
body interaction. The asymptotic convergence of the nuclear structure problem is not changed by 
including Jastrow type two-body correlation functions in the trial wavefunction [15 . The onset of 
asymptotic convergence occurs, however, at a much smaller value of an Af than for the case without 
correlation functions and the convergence is to the same final value of the lowest variational energy 
[3 HOI HI] . Therefore much less computational resources are required to get the answer. 

With the HO basis in the nuclear structure problem, convergence has been discussed, in practice, 
with an emphasis on obtaining those parameters which appear linearly in the trial function (i.e. 
convergence with Af) . Sometimes for each Af the non-linear parameter fku is varied to obtain the 
minimal energy |15j and then the convergence with Af is examined. Sometimes huj is simply fixed at 
a value which gives the fastest convergence in Af [5] . More recently, in the context of no-core shell 
model (NCSM) calculations and no-core full configuration (NCFC) calculations with smooth 'bare' 
potentials, one sees figures or tables in which one of the variational variables of (Af, Huj) is held 
fixed and the variational energy displayed with respect to the other. This practice is helpful for the 
following reason. Optimum values for the parameters that enter linearly can be obtained by solving 
a matrix eigenvalue problem. But the optimum value of the nonlinear parameter must in principle 
be obtained by, for example, numerical minimization which could be difficult as the algorithm could 
easily miss the global minimum and get trapped in a local minima. The plots one sees in the nuclear 
physics literature show that 1) for small bases a change in the non-linear parameter Huj can make a 
dramatic change in the variational estimate of the ground state energy and 2) the dependence on 
the nonlinear parameter decreases as the basis size increases. These observations seem to vitiate 
the need for an extensive numerical minimization by varying fku [32 . These observations have 
inspired definite (and differing) prescriptions for convergence and extrapolation. It is the purpose 
of this study to suggest that effective field theory (EFT) concepts of ultraviolet (uv) and infrared 
(ir) cutoffs provide an alternative useful way to think about convergence and a physically motivated 
prescription for extrapolation of (necessarily truncated) results in the model space (elucidated in 
Section 2) of the trial wavefunction to the full Hilbert space. 

The paper is organized as follows. In section 2 we briefly describe expansion schemes in HO 
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functions. This expansion technique still retains the variational character described above. We 
employ realistic smooth nucleon-nucleon potentials (JISP16 [33] and Idaho N 3 LO [34]) which have 
also been used by other authors without renormalization for A < 6 (Ref. [T2] and Ref. [TT] 
respectively). None of the discussion in section 2 is new, but it paves the way for section 3 in 
which we suggest a convergence analysis based upon the uv and ir cutoffs introduced in Ref. 35J 
in the context of an EFT framework. Section 4 is devoted to tests and examples of this new 
convergence scheme and section 5 contains a summary and outlook. 

II. EXPANSION IN A FINITE BASIS OF HARMONIC OSCILLATOR FUNCTIONS 

Here we indicate the workings of the finite HO basis calculations performed and refer the reader 
to a very useful review article [T3] on the no-core shell model (NCSM) for further details and 
references to the literature. In these no-core approaches, all the nucleons are considered active, so 
there is no inert core as in standard shell model calculations; hence the "no-core" in the name. TV TV 
potentials with strong short-range repulsions and the concomitant high-momentum components 
do not lend themselves well to a HO basis expansion, as was well appreciated fifty years ago [5]. 
A "renormalization" of the Hamiltonian is often made by constructing an effective interaction 
(dependent upon the basis cutoff Af and upon Huj) by means of a unitary transformation due to 
Lee and Suzuki |13j . This procedure generates effective many-body interactions that are often 
neglected [36 . This neglect destroys the variational nature of a NCSM calculation. We instead 
choose "soft" potentials (JISP16 [33J and Idaho N 3 LO [33]) which have also been used by other 
authors without renormalization for A < 6 (Ref. [H] and Ref. [IT] respectively), so that we 
can study convergence and extrapolation issues directly within a variational framework. NCSM 
calculations with these potentials are variational with the HO energy parameter hw and the basis 
truncation parameter Af as variational parameters [13] . Nomenclature has diverged somewhat since 
the advent of these smooth but still realistic potentials into a framework (NCSM) which originally 
included renormalization of the NN potential. Sometimes one reads about "NCSM calculations 
with unmodified or 'bare' potentials" [TTJ[T3], or "the no-core full configuration (NCFC) method" 
|12j . or simply "we use the basis of the no-core shell model (NCSM)" [37]. All these phrases refer 
to retaining the original interaction (without renormalization) within the model space. Nor do we 
renormalize the interaction in our study. 

To study these convergence issues we mostly employ the Idaho N 3 LO NN potential which is 
inspired by chiral perturbation theory and fits the two body data quite well [34] . It is composed 
of contact terms and irreducible pion-exchange expressions multiplied by a regulator function de- 
signed to smoothly cut off high-momentum components in accordance with the low-momentum 



7 



expansion idea of chiral perturbation theory. The version we use has the high-momentum cutoff 
of the regulator set at 500 MeV/c. The Idaho N 3 LO potential is a rather soft one, with heav- 
ily reduced high-momentum components as compared to earlier realistic NN potentials with a 
strongly repulsive core. Alternatively, in coordinate space, the Yukawa singularity at the origin is 
regulated away so that this potential would be considered "smooth" by Delves and Schneider and 
the convergence in J\f would be expected to be exponential [21 GO] • Even without the construction 
of an effective interaction, convergence with the Idaho N 3 LO NN potential is exponential, as nu- 
merous studies have shown [TTJ[37]. Nevertheless, it has been useful to simplify and reduce the 
high-momentum components of this and other phenomenological potentials further by means of 
the similarity renormalization group evolution |37j . Such a softening transformation is imperative 
for heaver nuclei (A > 6) and/or if three- nucleon forces are included in the Hamiltonian [3"?Tl39| . 

The second NN interaction we employ is JISP16 [33], a nonlocal separable potential whose 
form factors are HO wavef unctions. It is constructed by means of the J-matrix version of inverse 
scattering theory. The matrix of the NN potential in the oscillator basis is obtained for each 
partial wave independently, so the NN interaction is a set of potential matrices for different 
partial waves [40] . These matrices reproduce the experimental AW scattering data and properties 
of the deuteron to high precision. Once the inherent ambiguity of this method is eliminated by a 
plausible phenomenological ansatz, the scattering wavefunctions are very close to the ones provided 
by meson exchange "second-generation" NN potentials [41] . As for the name of this potential, 
JISP refers to J-matrix Inverse Scattering Potential and version "16" has had phase-equivalent 
unitary transformations applied to selected partial waves so that the resulting interaction continues 
to describe two-body data well. Selected partial waves are tuned to provide good descriptions of 3 H 
binding, the low-lying spectra of 6 Li and the binding energy of 16 [33] . The virtue of this potential 
is that it is also "soft". Although nonlocal and not really fitting into Delves classification, it is 
not surprising that variational calculations with this NN interaction also converge exponentially 
with J\f [12 since the HO form factors of this separable potential are gaussians multiplied by 
polynomials in the radial coordinate. (It is noteworthy that JISP16 in the HH basis also converges 
exponentially in K [42 , as would be expected by Schneider [10]). 

We use a HO basis that allows preservation of translational invariance of the nuclear self-bound 
system. Translational invariance is automatic if the radial HO wavefunction depends on relative, 
or Jacobi, coordinates as was done in Refs. [5, 7I-I10J . Antisymmetrization (or symmetrization for 
the alpha particle models of [TTllOj) of the basis is necessary and described in Refs. [T3] and [4"5] . 
Antisymmetrization in a Jacobi basis becomes analytically and computationally forbidding as the 
number of nucleons increases beyond four or five. For this reason these calculations are alternatively 
made with antisymmetrized wavefunctions constructed as Slater determinants of single-nucleon 
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wavefunctions depending on single-nucleon coordinates. This choice loses translational invariance 
since, in effect, one has defined a point in space from which all single-particle coordinates are 
defined. Translational invariance is restored by using the "Lawson method" [H] to be described 
shortly. The gain of this choice is that one can use technology developed and/or adapted for 
NCSM, such as the parallel-processor codes "Many-Fermion Dynamics — nuclear" (MFDn) 45j 
and the No-Core Shell Model Slater Determinant Code [46]. These codes set up the many-body 
basis space, evaluate the many-body Hamiltonian matrix, obtain the low-lying eigenvalues and 
eigenvectors using the Lanczos algorithm, and evaluate a suite of expectation values using the 
eigenvectors. 

The Slater determinant basis is often defined in the "m-scheme" where each HO single-particle 
state has its orbital and spin angular momenta coupled to good total angular momentum, j, and 
its magnetic projection, m. The many-body basis states for a given total number of nucleons A 
are Slater determinants in this HO basis and are limited by the imposed symmetries — parity, 
charge and total angular momentum projection (M), as well as by AT. In the natural parity cases 
for even nuclei M = 0, enables the simultaneous calculation of the entire low-lying spectrum for 
that parity and the chosen TV. 

The use of this specially constructed Slater determinant basis results in eigenstates of a trans- 
lationally invariant Hamiltonian (supplemented by a suitable constraint term) that factorize as 
products of a wavefunction depending on relative coordinates and a wavefunction depending on 
the CM coordinates. This is true for a particular truncation of the basis: a maximum of the sum of 
all HO excitations, i.e. X)iLi(2 n i + ^i) < N totmax , where n^,^ are the HO quantum numbers corre- 
sponding to the harmonic oscillators associated with the single-nucleon coordinates and N totma x is 
an example of the generic AT of the Introduction. Note that this truncation is on the level of total 
energy quanta ("total-energy-cut space"), which is different from the CI calculations used in atomic 
and molecular problems, which are often truncated on the single-particle level ("single-particle-cut 
space"). 

The precise method of achieving the factorization of the CM and intrinsic components of the 
many-body wavefunction follows a standard approach, sometimes referred to as the "Lawson 
method" [44] . In this method, one selects the many-body basis space in the manner described 
above with Af = N to tmax and adds a Lagrange multiplier term to the many-body Hamiltonian 
P{Hcm — ^hoj) where Hqm is the HO Hamiltonian for the CM motion. With j3 chosen positive 
(10 is a typical value), one separates the states of lowest CM motion (05i) from the states with 
excited CM motion by a scale of order f3huj. The resulting low- lying states have wavefunctions 
that then have the desired factorized form. We checked, for the two cases A = 3 and A = 4, that 
the codes manyeff (33] which use Jacobi coordinates and No-Core Shell Model Slater Determinant 
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Code [15] based upon single-nucleon coordinates gave the same eigenvalues for the same values 
of Af = Ntotmax and Hlo, indicating that the Lawson method is satisfactory for the calculations 
in single-particle coordinates. Some details of this check will be given in Section IV which gives 
results. 

Now we continue the discussion of the (total-energy-cut space) truncation parameter Af of the 
HO basis expansion of the many-body system. Usually, instead of truncating the sum of all 
HO excitations Af — N totmax > X^iLi(2 n i + h); one introduces the truncation parameter N max . 
N m ax is the maximum number of oscillator quanta shared by all nucleons above the lowest HO 
configuration allowed by the Pauli-exclusion principle for the chosen nucleus. We label the HO 

shells by energy quanta TV = (2n + I), where n = 0, 1, 2, . . . and I = 0, 1, 2, Thus, for example, 

in 6 He a truncation at N max = 4 would allow one neutron to occupy the TV — 5 HO shell, the 
other "valence" neutron would remain in the TV = 1 shell, and the remaining 4 nucleons remain 
in the (filled) first shell labeled by TV = 0. Alternatively, the two valence neutrons could occupy 
the TV — 3 shell and the remaining 4 nucleons stay in the TV — shell. In both cases (and for all 
other combinatorics) in 6 He N max = N to tmax — 2. Similarly for other p-shell nuclei one can work 
out that N to t m ax and N max differ, e.g. for 12 C, N max = N to tmax ™ 8, etc. However, for the s-shell 
nuclei 2 H, 3 H, 3 He, and 4 He N rnax = N tot max- 

Later on we will want to identify parameters of the model space (with the dimensions of mo- 
menta) which refer, not to the many-body system, but to the properties of the HO single-particle 
states. If the highest HO single-particle (SP) state of the lowest HO configuration allowed by the 
Pauli-exclusion principle has TV HO quanta, then N max + TV = N . Since N max is the maximum 
of the total HO quanta above the minimal HO configuration, we can have at most one nucleon in 
such a highest HO SP state with TV quanta. Note that N max characterizes the many-body basis 
space, whereas N is a label of the corresponding highest single-particle orbital. To find the value 
of the single-particle label N, we need to determine the highest occupied SP state in a given N max 
truncation. One gives all the available N max quanta to a single nucleon. Consider again a 6 He 
basis truncated at N max — 4; both valence neutrons occupy the Op (Nq = I) shell in the lowest 
energy many-body configuration. Assigning a single neutron the entire N max = 4 quanta means 
that, as before, the highest occupied SP state is in the N — 5 shell. On the other hand, the highest 
occupied orbital of the closed s-shell nucleus 4 He has TVo = so that TV = N max . 

III. CONVERGENCE IN UV AND IR VARIABLES 

We begin by thinking of the finite single-particle basis space defined by TV and hcj as a model 
space characterized by two momenta associated with the basis functions themselves. In the HO 
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basis, we define A = y // rriN(N + 3/2)hw as the momentum (in units of MeV/c) associated with the 
energy of the highest HO level. The nucleon mass is m^v = 938.92 MeV. To arrive at this definition 
one applies the virial theorem to this highest HO level to establish kinetic energy as one half the 
total energy (i.e., (N + 3/2)hui ) and solves the non-relativistic dispersion relation for A. This sets 
one of two cutoffs for the model space of a calculation. Energy, momentum and length scales are 
related, according to Heisenberg's uncertainty principle. The higher the energy or momentum scale 
we may reach, the lower the length scale we may probe. Thus, the usual definition of an ultraviolet 
cutoff A in the continuum has been extended to discrete HO states. It is then quite natural to 
interpret the behavior of the variational energy of the system with addition of more basis states as 
the behavior of this observable with the variation of the ultraviolet cutoff A. Above a certain value 
of A one expects this running of the observable with A to "start to behave" so that this behavior 
can be used to extrapolate to the exact answer. However, the model space has another scale which 
motivates a second cutoff; the energy scale of hui itself. Because the energy levels of a particle in 
a HO potential are quantized in units of fuv, the minimum allowed momentum difference between 
single-particle orbitals is A = y/m^huj and that has been taken to be an infrared cutoff |35j . That 
is, there is a low-momentum cutoff A = h/b corresponding to the minimal accessible non-zero 
momentum (here b = ^/ m j^ plays the role of a characteristic length of the HO potential and basis 
functions). Note however that there is no external confining HO potential in place. Instead the 
only huj dependence is due to the scale parameter of the underlying HO basis. 

The energies of a many-body system in the truncated model space will differ from those calculated 
as the basis size increases without limit (JV — > oo). This is because the system is in effect confined 
within a finite (coordinate space) volume characterized by the finite value of b intrinsic to the HO 
basis. The "walls" of the volume confining the interacting system spread apart and the volume 
increases to the infinite limit as A -> and b — > oo. Thus it is as necessary to extrapolate the 
low momentum results obtained with a truncated basis with a given b or fvjj as it is to ensure 
that the ultraviolet cutoff is high enough for a converged result. These energy level shifts in a 
large enclosure have long been studied [47]; most recently with the explicit EFT calculation of 
a triton in a cubic box allowing the edge lengths to become large (and the associated ir cutoff 
due to momentum quantization in the box going towards zero) |48) . There it was shown that as 
long as the infrared cutoff was small compared to the ultraviolet momentum cutoff appearing in 
the "pionless" EFT, the ultraviolet behavior of the triton amplitudes was unaffected by the finite 
volume. More importantly, from our point of view of desiring extrapolation guidance, this result 
means that calculations in a finite volume can confidently be applied to the infinite volume (or 
complete model space) limit. Similar conclusions can be drawn from the ongoing studies of systems 
of two and three nucleons trapped in a HO potential with interactions from pionless EFT and this 
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definition (A = ^/rn^huj) of the infrared cutoff [39] . 

Other studies define the ir cutoff as the infrared momentum which corresponds to the maximal 
radial extent needed to encompass the many-body system we are attempting to describe by the 
finite basis space (or model space). These studies find it natural to define the ir cutoff by A sc = 

{wiNhuj) / {N + 3/2) j37J[SD]. Note that A sc is the inverse of the root-mean-square (rms) radius 
of the highest single-particle state in the basis; {r 2 ) 1 / 2 = by/ N + 3/2. We distinguish the two 
definitions by denoting the first (historically) definition by A and the second definition by A sc 
because of its scaling properties demonstrated in the next Section. This latter ir variable A sc 
clearly goes to zero either i) as huj goes to zero at fixed TV or ii) as N becomes large for fixed ttw. 
It is the second limit which corresponds to the convergence theorems of the Introduction. In this 
latter limit X sc — > while A does not. 

The extension of the continuum ultraviolet cutoff to the discrete (and truncated) HO basis with 
the definition A = \JmN(N + 3/2)/kj seems unexceptional. An equally plausible alternative uv 
cutoff differs from A by only a scale change [5T], in striking contrast to the alternate definitions 
of the ir cutoff which have different functional forms. It is a goal of this work to determine 
the usefulness of the two rival definitions- A and A sc -of the infrared regulator (let us call it Xi r ) 
of the model space with parameters (N, Huj). From the beginning, it is clear that increasing 
A by increasing Hut in a fixed- N model space is not sufficient; doing so increases both of the 
putative infrared cutoffs as well because A = X^/N + 3/2 = X SC (N + 3/2) and one continues to 
effectively calculate in an effective confining volume which is getting smaller rather than larger. 
This confining volume is certainly removed by letting N — > oo, at fixed fku, because HO functions 
form a basis of the complete space. But increasing N without limit is computationally prohibitive, 
and furthermore will not shed much light on the question of whether in practical calculations one 
must take the ir cutoff to zero by taking huj — > (Aj r = A definition) or whether it is sufficient to 
allow hu be some larger value, perhaps near that used in traditional shell-model calculations, and 
let an increasing N take Xi r to small values, as it does with the definition Xi r = A sc . 
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FIG. 1: (Color online) Schematic view of a finite model space (limited by the basis truncation parameter 
N as described in the text), in which the uv and ir momentum cutoffs are arbitrary. To reach the full 
many-body Hilbert space, symbolized by the complete oval, one needs to let the uv cutoff — > oo and the ir 
cutoff — > 0. 

We are interested in the limit of large A and small X ir (see Figure 1). If one can establish that 
the cutoff dependences of the model space decrease with increasing A and decreasing \ r then 
one can i) remove the influence of the ir cutoff by extrapolating to the infrared limit for selected 
uv-cutoff values chosen to be above the uv nature of the potential and ii) if needed, extrapolate to 
the uv limit for selected ir cutoff values chosen by the size of the system modeled. We will show 
that such a program is possible. 

IV. RESULTS AND DISCUSSION 

We first display in a series of figures the running of the ground-state eigenvalue of a single 
nucleus, 3 H, on the truncated HO basis by holding one cutoff of (A, Aj r ) fixed and letting the 
other vary. Then we show that the trends noted hold for other light nuclei within the range of our 



13 




FIG. 2: (Color online) Dependence of the ground-state energy of 3 H (compared to a converged value-see 
text) upon the uv momentum cutoff A for different fixed A. The curves are fit to the calculated points. 

computer resources. Finally we discuss extrapolation procedures. 

These 3 H calculations were made for N < 36 and values of hut as appropriate for the chosen cutoff 
value. For N > 16, we used the code manyeff [13] which uses Jacobi coordinates and the No-Core 
Shell Model Slater Determinant Code [IS] which use single-particle coordinates for smaller N. We 
checked that the codes gave the same eigenvalues for overlapping values of N, indicating that the 
Lawson method satisfactorily restores translational invariance to ground-state energy calculations 
in single-particle coordinates. For example, the ground state energy of 3 H with the Idaho N 3 LO 
NN potential at N max = 16 and huj = 49.2968 MeV is (-7.3378, -7.3385) MeV for the (Jacobi, 
single-particle) basis choice. 

In Figure 2 and the following figures, \AE/E\ is defined as \(E(A,X ir ) — E)/E\ where E 
reflects a consensus ground-state energy from benchmark calculations with this NN potential, 
this nucleus, and different few-body methods. The accepted value for the ground state of 
3 H with this potential is —7.855 MeV from a 34 channel Faddeev calculation 34J, —7.854 
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MeV from a hyperspherical harmonics expansion [S2], and —7.85(1) from a NCSM calculation 
[IT]. All \AE/E\, starting with Figure 2, will follow some trajectory (trajectory's shape not 
predicted). For the choice of Figure 2, X ir = A = ^/rrv^huj, \AE/E\ — > as A increases. 
Fixed hui implies N alone increases to drive A — > oo, \ sc —¥ simultaneously. The linear fit 
on a semi-log plot is extracted from the data. This fit implies \AE/E\ ~ B cxp(— A/A re /(A)), 
where B is approximately constant and cA re /(A) ~ 30fku for Hui > 45 MeV. Note this is 
■A-re/M an( i n °t A re j(A sc ), i.e. with fixed A, A re f is a constant. On the other hand, for fixed A, 
a smaller A implies a smaller \ AE/E\ since more of the infrared region is included in the calculation. 
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FIG. 3: (Color online) Dependence of the ground-state energy of 3 H (compared to a converged value-see 
text) upon the ir momentum cutoff A for fixed A. 

In Figure 3 we hold fixed the uv cutoff of (A, A^ r ) to display the running of \AE/E\ upon the 
suggested ir cutoff A. For fixed A, a larger A implies a smaller | AE/E\ since more of the uv region 
is included in the calculation. But we immediately see a qualitative change in the curves between 
the transition A = 700 MeV and A = 800 MeV; for smaller A, \AE/E\ does not go to zero as the 
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ir cutoff is lowered and more of the infrared region is included in the calculation. This behavior 
suggests that \AE/E\ does not go to zero unless A > A NN , where A^* is some uv regulator scale 
of the NN interaction itself. From this figure one estimates A*" ~ 800 MeV/c for the Idaho 
N 3 LO interaction. 

Yet the description of this interaction in the literature says that the version we use has the 
high-momentum cutoff of the regulator set at An^lo — 500 MeV/c [34 . This does not mean 
that the interaction has a sharp cutoff at exactly 500 MeV/c, since the terms in the Idaho N 3 LO 
interaction are actually regulated by an exponentially suppressed term of the form 



In this expression, p and p' denote the magnitude of the initial and final nucleon momenta of this 
non-local potential in the center-of-mass frame and n > 2. Because the cutoff is not sharp, it 
should not be surprising that one has not exhausted the uv physics of this interaction for values of 
single-particle A somewhat greater than 500 MeV/c. Note that this form of the regulator allows 
momentum transfers (p — p 1 ) to achieve values in the range up to 2A.N3LO- Can one make an 
estimate of the uv regulator scale of the Idaho N 3 LO interaction which is more appropriate to the 
discrete HO basis of this study? An emulation of this interaction in a harmonic oscillator basis 
uses Huj — 30 MeV and N max = N = 20 [23]. Nucleon- nucleon interactions are defined in the 
relative coordinates of the two-body system so one should calculate A 7 ™ = y/m(N + 3/2)hw with 
the reduced mass m rather than the nucleon mass appropriate for the single-particle states 
of the model space. Taking this factor into account, the successful emulation of the Idaho N 3 LO 
interaction in a HO basis suggests that A NN ~ 780 MeV/c, consistent with the figure. 

For A < A NN there will be missing contributions of size |(A — A NN )/A NN \ so "plateaus" develop 
as A — > revealing this missing contribution to \AE/E\. The "plateau" is not flat as A —> and, 
indeed, rises significantly with decreasing A < A NN . This suggests that corrections are needed to 
A and A which are presently defined only to leading order in A/ A. We hope to learn if higher-order 
corrections can be determined by the data in a future study. 
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FIG. 4: (Color online) Dependence of the ground-state energy of 3 H (compared to a converged value-see 
text) upon the uv momentum cutoff A for different values of the ir momentum cutoff A sc . Curves are not 
fits but simple point-to-point line segments to guide the eye. 

Now we turn to the second pair of cutoffs of (A, A, r ) and display in Figure 4 the analogue of 
Figure 2 except that this time \i r = X sc = yjm^hjjj (N + 3/2). For fixed X sc , / E\ does 
not go to zero with increasing A, and indeed even appears to rise for fixed X sc > 35 MeV/c 
and A > 800 MeV/c. Such a plateau-like behavior was attributed in Figure 3 to a uv regulator 
scale characteristic of the NN interaction. Can the behavior of Figure 4 also be explained by a 
"missing contributions" argument; i.e. an argument based upon A sc < A^,^ where A^^ is a second 
characteristic ir regulator scale implicit in the NN interaction itself? One can envisage such an ir 
cutoff as related to the lowest energy configuration that the NN potential could be expected to 
describe. For example, the inverse of the np triplet scattering length of 5.42 fm corresponds to a 
low-energy cutoff of about 36 MeV/c. Realistic NN potentials such as Idaho N 3 LO and JISP16 
do fit these low-energy scattering parameters well. The previously mentioned emulation of the 
Idaho N 3 LO interaction in a harmonic oscillator basis [53] has X^ C N ~ 36 MeV/c. As we shall 
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see later, the fit to low-energy NN data of JISP16 implies a <~ 63 MeV/c. The factor of 

two difference may simply be a reflection of the fitting procedures and appears to be within the 
range of uncertainty of our argument. In any event, the behavior of our results in Figure 4 is not 
inconsistent with this concept of an inherent ir regulator scale implicit in the NN interaction. 

Having introduced a scale X^ C N , we continue our discussion of Figure 4. As fixed A sc requires 
hui/N to be constant and N < 36, small values of fixed X sc are linked with small values of A. Having 
said that, we see that \AE/E\ — > with increasing A for X sc — 20 McV/c and perhaps A sc = 28 
McV/c for the values of A available to the calculation. At fixed X sc > Xf c N and increasing A, once 
A > A^^, a "plateau" will develop since no new contributions to \AE/E\ exist for A > A NN . 
In the figure, the plateau appears to start at A ~ 700 — 900 MeV/c, consistent with the estimate 
of A NN ~ 780 MeV/c for this NN interaction. The plateau in \AE/E\ for larger fixed A sc is 
higher than the plateau for small fixed X sc since more contributions to \AE/E\ are missing from 
the infrared region. Again we observe the plateau rises with increasing A and this behavior may be 
a sign that corrections are needed to A and A sc which are presently defined only to leading order 
in A sc /A. However, for A sc < X^ C N and A < A NN the results converge to a single curve at the 
left of this figure. It is remarkable that this curve persists to quite low A values. This means that 
\AE/E\ becomes insensitive to X sc for low A if X sc is low enough. Later on we will demonstrate 
that this curve can be quite well described by a Gaussian, a result which persists for other s-shcll 
nuclei. But we will see in the next figure that one has not yet captured the uv region at these low 
values of A. 

Figure 5 is the analogue to Figure 3: only the variable on the x-axis changes from A to A sc = A 2 /A. 
For A < A NN ~ 780 MeV/c the missing contributions and resulting "plateaus" are as evident as 
in Figure 3. The tendency of these plateaus to rise as A sc — > again suggests a refinement is 
needed to this first-order definition of the cutoffs. Around A <~ 600 MeV/c and above the plot of 
\AE/E\ versus A sc in Figure 5 begins to suggest a universal pattern, especially at large A sc . For 
A <~ 800 MeV/c and above the pattern defines a universal curve for all values of A sc . This is the 
region where A > A^^indicating that nearly all of the ultraviolet physics set by the potential has 
been captured. Such a universal curve suggests that A sc could be used for extrapolation to the ir 
limit, provided that A is kept large enough to capture the uv region of the calculation. Figure 5 
is also the motivation for our appellation A sc , which we read as "lambda scaling", since this figure 
exhibits the attractive scaling properties of this regulator. 

For Figure 6, we take advantage of the "saturation" of the uv region by binning all results with 
A > 800 MeV/c. They do indeed fall on a universal curve for each nucleus shown, indicating that 
one can use this universal behavior for an extrapolation which is independent of A. The curves 
are offset for the three nuclei but otherwise appear similar. The points can be fit by the function 
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FIG. 5: (Color online) Dependence of the ground-state energy of 3 H (compared to a converged value-see 
text) upon the ir momentum cutoff A sc for fixed A. 

y = acxp(-b/\ sc ) with b « 20 - 40 MeV/c at A=2,3,4. 

Finally we utilize the scaling behavior displayed on Figures 5 and 6 to suggest an extrapolation 
procedure which we demonstrate in Figure 7, again concentrating on 3 H and the Idaho N 3 LO 
potential. The extrapolation is performed by a fit of an exponential plus a constant to each set of 
results at fixed A. That is, we fit the ground state energy with three adjustable parameters using 
the relation E gs (X sc ) = aexp(— b/X sc ) + E gs (\ sc = 0). It should be noted that our five extrapola- 
tions in Figure 7 employ an exponential function whose argument 1/A SC = \J (N + 3/2)/(mjv^w) 
is proportional to \J N /(hui) and is therefore distinct from the popular extrapolation with an expo- 
nential in N max (= N for this s-shell case) [TTHT31 [571 138 , 53, 54]. The mean and standard deviation 
of the five values of E gs (X sc = 0) were —7.8511 MeV and 0.0011 MeV, respectively, as suggested 
by Figure 7 in which the overlap of the five separate curves cannot be discerned. For calibration 
of our scheme, we recall that the accepted value for the ground state of 3 H with this potential 
is —7.855 MeV from a 34 channel Faddeev calculation |34| . —7.854 MeV from a hyperspherical 
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FIG. 6: (Color online) Dependence of the ground-state energy of three s-shell nuclei (compared to a 
converged value-see text) upon the ir momentum cutoff A sc for A above the w 780 MeV/c set by the 

NN potential . 

harmonics expansion [55], and —7.85(1) from a NCSM calculation |11) . 

This extrapolation procedure with a large and fixed A = \J m^iN + 3/2)huj and taking A sc = 
^/mjyhui/(N + 3/2) toward the smallest value allowed by computational limitations treats both TV 
and huj on an equal basis. For example, the extrapolation at fixed A = 1200 MeV/c employs values 
of Hlu from 41 to 65 MeV and N = 22 - 36. The one at fixed A = 800 MeV/c employs values of 
huj from 18 to 44 MeV and N = 14 — 36. The curves of Figure 7 encompass values of A sc between 
20 and 52 MeV/c. We attempted to quantify the spread in extrapolated values by fitting only 
segments of the curves of this figure. Recall that the smallest value of A sc requires the largest N. 
Fits to the segment from A sc = 20 MeV/c to X sc = 40 MeV/c (always for the five displayed values 
of fixed A) resulted in a mean of —7.8523 MeV and standard deviation of 0.0008 MeV. Cutting out 
the left hand parts of the curves and fitting only from A sc = 30 MeV/c to A sc = 55 MeV/c gave a 
mean of —7.8498 MeV and standard deviation of 0.0022 MeV. For both these trials a rather large 
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FIG. 7: (Color online) The ground state energy of 3 H calculated at five fixed values of A and variable X sc . 
The curves are fits to the points and the functions fitted are used to extrapolate to the ir limit X ac — 0. 

N was needed, ranging from 14 to 36 but the extrapolation is quite stable. In contrast, values 
of A sc higher than those shown in Figure 7, namely from A sc = 50 MeV/c to A sc = 85 MeV/c, 
require fewer computational resources (N — 8 — 22). The extrapolations have a mean and standard 
deviation of —7.792 MeV and 0.042 MeV, still not so far away from the accepted value of —7.85 
MeV. 

In Figure 8 we return to one of the curves of Figure 4 to examine the dependence of the binding 
energy of three nuclei upon the uv regulator when the putative ir regulator is held fixed (A sc = 55 
MeV/c). At fixed X sc > X^ C N ~ 36 MeV/c, and increasing A, once A > A NN , a "plateau" will 
develop since no new contributions to \AE/E\ exist for A > A"" ~ 780 MeV/c. The new feature 
of this figure is that the "plateau" of the nucleus 2 H is above that of 3 H (taken from figure 4) which 
is in turn above that of 4 He. This suggests that Ajviv is not the only regulator scale needed to 
explain the dependencies upon A and A sc . Figure 8 introduces another scale - the role of binding 
momentum (Q) of a nucleus. The scale Q has been used recently in EFT treatments of pion- 
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FIG. 8: (Color online) Dependence of the ground-state energy of three s-shell nuclei (compared to a 
converged value-see text) upon the uv momentum cutoff A for A sc above the X^ C N « 36 MeV/c set by the 
NN potential. Curves are not fits but spline interpolations to guide the eye. 

deuteron scattering at threshold [55] . The idea is to take the small binding energy of the deuteron 
explicitly into account as one attempts to develop a consistent power counting for an EFT of pion- 
nucleus scattering lengths. The extension of the definition of Q to more massive nuclei can take 
alternate forms: Q — ^2m^{E / A) where E/A is the binding energy per nucleon, or Q — y/2/ie 
where \x is the reduced mass of a single nucleon with respect to the rest of the nucleons in the 
nucleus and e is the binding energy with respect to the first breakup channel |56j . Clearly the two 
definitions coincide for the deuteron and for the light nuclei considered here both definitions give 
similar estimates. For definitiveness, we calculate Q according to the formula Q = \/2/j,e. This 
calculation gives Q( 2 H) = 46 MeV/c, Q( 3 H)= 88 MeV/c, Q( 3 He)= 83 MeV/c , Q( 4 He)= 167 
MeV/c, and Q( 6 He) = 39 MeV/c. The binding momentum of 6 He is comparable to that of the 
deuteron because the first breakup channel into 4 He+2n is only about 1 MeV above the ground 
state. 
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The fractional error plotted in Figure 8 appears to rise slightly from a minimum at A > A^" ~ 
780 MeV/c as the uv cutoff A increases for each of the three nuclei. For example in the 2 H 
calculation the X sc cutoff relative to the deuteron binding momentum is X sc /Q = 1.2 and the 
error is rather high, rising from a minimum of about 25%. The triton is more bound so that the 
ratio is X sc /Q = 0.62 and the minimum error is 9%. The calculation of the tightly bound 4 He 
(^sc/Q = 0.33) has the smallest error of less than 2%, but even that error appears to rise as A 
increases to the limits of the present calculation. That is, as Q increases at fixed A sc (at high 
enough A) the error due to the A sc cutoff is lower. It is natural to expect that the many-body 
dynamics enters at some level and sets additional scales beyond the TVA^-interaction scales. 
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FIG. 9: (Color online) Dependence of the ground-state energy of three s-shell nuclei (compared to a 
converged value-see text) upon the uv momentum cutoff A for A sc below the A^' « 36 MeV/c set by the 
NN potential. The data are fit to Gaussians. 

A final comparison of three s-shell nuclei selects the low end of the ir region where X sc is at or 
below the ir cutoff suggested by the potential. In Figure 9 all momenta are scaled by the binding 
momentum Q of the considered nucleus in order to put them on the same plot. For such low 
momenta A sc , \AE/E\ does go to zero with increasing A because X sc < Xf c N , where is the 
second regulator scale of the NN interaction itself. For 3 H, \ SC /Q ~ 0.23 corresponds to A sc = 20 
MeV/c; the curve can be directly compared with the analogous curve (black online) in Figure 4. 
For 2 H A sc = 10 MeV/c and for 4 He X sc — 40 MeV/c, all values of X sc are near or below the second 
(ir) regulator scale of the Idaho N3LO potential suggested to be ~ 36 MeV/c. The largest value 
of A plotted is - 861 MeV/c for 2 H, - 746 MeV/c for 3 H, and - 780 MeV/c for 4 He, all values 
of A are near or just above the uv regulator of the Idaho N3LO potential suggested to be ~ 780 
MeV/c. So this is a plot of low A for all the nuclei portrayed. The "high" A tails of these curves 
can be fit by Gaussians (shifted from the origin) in the variable A/Q. This figure would seem to 
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nicely illustrate the expectations of the theorems of the 1970's [5J |5D] that the asymptotic rate of 
convergence does not depend upon the number of particles. 
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FIG. 10: (color online) The ground state energy of 6 He calculated at five approximate values of A (see 
text) and integer values of hu). The curves are fits to the points and the functions fitted are used to 
extrapolate to the ir limit A = as in Figure 7. 

In Figures 4 through 9 we have displayed the features of our results as functions of the pair 
of cutoffs of (A, Ai r ) where A; r = A sc = \J m^tuo / (N + 3/2), and demonstrated an extrapolation 
procedure to the uv and ir limits. Yet Figure 3 suggests that an extrapolation to the infrared limit 
could equally well be made by taking A — > for a fixed large A. In Figure 10, we demonstrate the 
features of such an extrapolation by using published results of the JISP16 interaction for the halo 
nucleus 6 He [12] ■ fts binding momentum Q( 6 He) = 39 MeV/c is comparable to that of the loosely 
bound deuteron and should provide a severe test of any infrared extrapolation. 

A second reason for considering 6 He is that it has been studied extensively with the JISP16 NN 
interaction, both with HO expansion and HH expansion techniques |57j . The same set of ground 
state energy eigenvalues as those plotted in Figure 10 yields an extrapolated value of —28.76(9) MeV 



25 



using "extrapolation A" of Ref. [HI [55] . In "extrapolation A" one lets the variable N — > oo with a 
selection procedure for values of Huj as explained in Ref. [12] . If one refers to the published results 
[T!2] where the largest N is 15 one finds an "extrapolation B" method which lets N — > oo at fixed tkJ 
to obtain —28.69(5) MeV in good agreement with the value of — 28.68(12)MeV from "extrapolation 
A". A hyperspherical harmonics expansion calculation of 6 He with the JISP16 potential finds an 
extrapolated value of -28.70(13) MeV [42] . This is increased by about 200-300 keV to -28.96(3) 
MeV by a "hyperspherical harmonics effective interaction" technique which requires fewer terms 
in K to reach asymptotic convergence but loses variational character because the induced many- 
particle interactions are dropped from the effective interaction |23) . 

The results of Ref. |12j were obtained with an antisymmetrized many-body wavefunction con- 
structed as as sum of Slater determinants of single-nucleon wavefunctions depending on singlc- 
nucleon coordinates (and the Lawson method to isolate CM effects) on a mesh of integer (N, fuS). 
The value of N is, by definition an integer and values of the non-linear variational parameter Huj 
were chosen to be increments of 2.5 MeV between 10 and 40 MeV. To show that the familiar integer 
values of Hoj from Ref. |12| could be directly used in the extrapolation procedures suggested here 
we mapped the ground state energy eigenvalues onto the variables (A, tiw = A 2 /mjv) rather than 
onto the variables (A, y/huj = A/y^mjv). The largest value of N was 17 (N max = 16 for this p-shell 
nucleus) [58) . The extrapolation of Figure 10 is performed by a fit of an exponential plus a constant 
to the set of results at fixed A. The resulting A's are not then strictly fixed but each point plotted 
corresponds to a value of A constant to within 2 — 5% of the central value indicated on the graph. 
The important S'-wave parts of the JISP16 potential are fit to the data in a space of N = 8 and 
hco = 40 MeV. Therefore this potential has NN regulator scales of X^ C N ~ 63 MeV/c and A NN ~ 
600 MeV/c. (The value of X NN associated with this potential is about 200 MeV/c, as can readily 
estimated from the legend of Figure 2). But JISP16 seems to be so soft that the ultraviolet region 
is already captured with A > 500 MeV/c, as shown by the top two curves of Figure 10. 

We fit the ground state energy with three adjustable parameters using the relation E gs {tno) — 
a exp(— 6 / ftw) + E gs (huj — 0) five times, once for each "fixed" value of A. It is readily seen that 
one can indeed make an ir extrapolation by sending huj — > with fixed A as first advocated in Ref. 
35j and that the five ir extrapolations are consistent. The spread in the five extrapolated values is 
about 500 keV or about 2% about the mean of —28.78 MeV. The standard deviation is 200 keV. 

A second (single) extrapolation of the 6 He data with A sc — > which uses all calculated energies 
where A > 510 MeV/c is shown in Figure 11. As in Figure 7, we fit the ground state energy with 
three adjustable parameters using the relation E gs (X sc ) = aexp(— b/X sc ) + E gs (X sc — 0). The 
extrapolated value is —28.68 MeV, which agrees well with Figure 10 and the other extrapolated 
results. The extrapolation prescription used in Figure 11 employs values of Huj from 15 to 40 MeV 
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and a range of N from 7-17. That is, all of the information (at A > 510 MeV/c) available from 
these calculations is used in the A sc extrapolation. How can one estimate an uncertainty from such 
a single extrapolation? Looking at the scatter of the points about the fitted curve is instructive but 
not quantitative. If we bin the 48 points of Figure 11 into the same bins of "constant" A as in Figure 
10, we find (not shown) five extrapolations with a mean of —28.58 MeV and standard deviation of 
0.06 MeV. Another possible way of breaking up this single extrapolation is more in the spirit of 
the earlier extrapolations of Delves and successors. If we map from (N, tuo) onto (A, A sc ) holding 
N fixed we get another set of extrapolations, those for N = 7,9, 11, 13, 15, 17 with concomitant 
smallest A sc = 60,50,41,36,32,28 MeV/c. Of these six extrapolations, only those with N > 13 
are consistent with the extrapolation which uses the full 41 points. This is to be expected, as one 
needs a large N before the convergence "starts to behave" . The mean of the three extrapolations 
with N > 13 is —28.54 MeV and standard deviation is 0.11 MeV. Concentrating only on large N 
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in this naive manner gives a worse extrapolation compared to accepted extrapolated ground state 
energies. It would appear that it is advantageous to take advantage of the scaling properties of 
\ sc for all values of the uv regulator large enough to capture the uv limit. In that case, as seen in 
Figure 11, even results with low N (and therefore large A sc ) can usefully stabilize and bound an 
extrapolation to the ir limit. A rough estimate of the uncertainties of this extrapolation of figure 
11 would then be -28.68(22) MeV. 

In conclusion, our extrapolations in the ir cutoff A of —28.78(50) MeV or the ir cutoff A sc of 
28.68(22) MeV are consistent with each other and with the independent calculations. 

V. SUMMARY AND OUTLOOK 

We reviewed the functional analysis theorems which describe variational calculations of many- 
body systems made with a trial function expanded in a complete set of known functions. According 
to these theorems the convergence properties of such a calculation are determined by the interaction 
and by the dimensionlcss number Af which determines the truncation at a finite number of basis 
functions. Among basis sets, harmonic oscillator (HO) functions are distinguished by ease of 
separation of relative and center of mass coordinates and by the dimensional parameter Huj which 
sets an intrinsic scale. Motivated by effective field theory studies, one can define quantities from 
(Af,Huj) forming ultraviolet (uv) and infrared(ir) momenta that act as cutoffs that characterize the 
model space just as does (Af, Huj). Extending both the uv cutoff to infinity and the ir cutoff to zero 
is prescribed for a converged calculation. There have been two alternate definitions of the ir cutoff; 
A = y/miqHuj and A sc = y/ra^Huj/ (N + 3/2). Note that A sc = A 2 /A where A = y/m,N(N + 2>/2)Huj 
is the uv cutoff as usually defined. We calculated the ground state energy of light nuclei with the 
"bare" and "soft" NN interactions Idaho N 3 LO and JISP16. We investigated the behaviors of the 
uv and ir regulators of model spaces used to describe 2 H, 3 H, 4 He and 6 He. 

We obtained fully converged eigenvalues for 2 H and 3 H which were in agreement with other 
calculations including those (e.g. Faddeev approach) obtained from a direct finite difference solu- 
tion of partial differential equations in many dimensions. These results could be used to examine 
the cutoff dependences of the model spaces (A, A) or (A, A sc ) as one cutoff was held fixed and 
the other approached its limit. The examination was based upon the ratio \AE/E\, defined as 
\(E(A, \ ir ) —E)/E\ where E is the fully converged ground state energy. Both pairs of cutoffs acted 
as expected when A was held fixed and X ir tended toward zero; \AE/E\ — > in Figures 3 and 5 
provided A exceeds a threshold set by the potential. On the other hand, in both figures drawing 
an imaginary vertical line at a fixed \ r which crosses the curves shows that the calculation gets 
better as A increases. It is in Figures 2 and 4 that the difference between the two versions of the 
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Ai r cutoff become evident. For the pair (A, A) in Figure 2 \AE/E\ — > as the uv cutoff increases 
for all values of A investigated. But in Figure 4 \AE/E\ actually rises as A — > oo if X sc is larger 
than another threshold value evidently set by the potential. As both cutoffs should be sent to their 
respective limits for a converged calculation, this behavior does not invalidate the identification of 
A sc with Ai r , but it does seem a little peculiar. Perhaps this behavior signals a need for higher 
order terms in A/ A in the definition of A sc 

In any event, we have introduced a practical extrapolation procedure with A — ¥ oo and \ r — > 
which can be used when the size of the basis exceeds the capacity of the computer resources as 
it does for He and 6 He and certainly will for any more massive nuclei. We have established that 
A does not need to be extrapolated to oo but if A > K NN set by the potential one can make the 
second extrapolation to zero with either ir cutoff A sc (see Figures 7 and 11) or A (see Figure 10). 
The choice of the scaling cutoff A sc is especially attractive as A need not be held constant but any 
A large enough can be used in the ir extrapolation. Unlike other extrapolation procedures the ones 
advocated in this paper treat the variational parameters Af and hu on an equal footing to extract 
the information available from sequences of calculations with model spaces described by (Af , Huj) 

For the future, we can envisage extending this extrapolation technique to calculating other prop- 
erties of nuclei, properties which may or may not be as amenable as are energy eigenvalues to the 
uv and ir regulators. The mis point matter radii and the Gamow- Teller matrix element (relevant 
to /3 decay) of light nuclei are important quantities to calculate reliably for these (and more massive 
nuclei) [42, 59]. In the nuclear structure folklore, r 2 and D z (the z component of the electric-dipole 
operator) are of long range and the full GT matrix element, including meson-exchange currents, is 
of medium range. The electric dipole polarizabilities of light nuclei are necessary in order to obtain 
accurate nuclear-polarization corrections for precisely measured transitions involving «S-waves in 
one-and two-electron atoms. The defining relation for the polarizability can be converted into a 
procedure which needs only bound-state quantities and involves the long-range dipole operator D 
[50] , A convergence analysis of a HO expansion, which lets N — > oo at fixed hui, for electric dipole 
polarizabilities of 3 H, 3 He and 4 He obtained faster convergence for lower hui than for the binding 
energy itself |61j . It would be interesting to learn how the procedure advocated here would work 
for these problems. These latter problems, often require not only converged ground state energies, 
but energies which agree with experiment. For that, NNN interactions arc considered necessary 
[62j , thereby leading to a need for more studies of the convergence and extrapolation concepts of 
this paper. 
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